COULOMB INTERFERENCE IN HIGH-ENERGY pp 
AND pp SCATTERING. 



V. A. PETROV a '[], E. PREDAZZI 6 -g and A. PROKUDIN^'Q 

(a) Institute For High Energy Physics, 
142281 Protvino, RUSSIA 

(b) Dipartimento di Fisica Teorica, 
Universitd Degli Studi Di Torino, 
Via Pietro Giuria 1, 10125 Torino, 
ITALY 
and 

Sezione INFN di Torino, 
ITALY 



An analysis of the Coulombic amplitude and its interference with the nuclear 
amplitude which is driven by the three-component Pomeron [Q is presented. It 
is shown that different approaches towards the Coulomb phase evaluation give 
approximately uniform result at all energies and the difference is negligible at 
RHIC and LHC energies. We show that the use of the amplitude which was 
fitted to accomodate nucleon data only (in the region 0.01 < \t\ < 14.5 (GeV 2 )) 
combined with the Coulomb amplitude, reproduces the existing data in the 
Coulomb interference domain quite accurately without any adjustement of the 
parameters. As a consequence, we predict the differential cross section in the 
region of the Coulomb nucleon interference for both RHIC and LHC energies. 



1 INTRODUCTION 



In Ref. [|X|] an eikonal model of a three-component Pomeron has been sug- 
gested and successfully used for describing the high energy pp nad pp data 
in the region of large momentum transfer 0.01 < \t\ < 14.5 (GeV 2 ) . In 
this paper we apply the model to the region of small momentum transfer 
< \t\ < 0.01 (GeV 2 ). 
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The problem is a proper account of the Coulomb interaction which is most 
important at the smallest \t\. The standard way to do this is to represent the 
whole scattering amplitude T(s, t) which is dominated by the Coulomb force 
at low momentum transfer and by the hadronic force at higher momentum 
transfer as 

T(s,t) = T N (s,t)+e^T c (s,t), (1) 
where if we normalize the scattering amplitude so that 

da _ \T(s,t)\ 2 

dt 167TS 2 ' 

the Born Coulomb amplitude for pp and pp scattering is 

T c (s } t)=T—. (3) 

The upper (lower) sign corresponds to the scattering of particles with the 
same (opposite) charges. T N (s, t) stands for purely strong interaction ampli- 
tude, and the phase <3> depends generally on energy, the momentum trasfer 
and on the properties of T . The study of the Coulomb nuclear interfer- 
ence is very important for extracting the real part of the strong interaction 
amplitude. 

The issue of Coulombic amplitude and its interference with the nuclear 
component was addressed in many papers in the past. West and Yennie || 
examined the Coulomb-nuclear interference using Feynman diagrams. R. 
Cahn H considered the same task in the framework of the Eikonal model; 
his results were quite convincing (though the modifications were very small 
compared to @) and proved that the Eikonal model is a very convenient 
basis for analyzing the Coulomb-nuclear interference. 

The Coulombic phase has atracted the attention of many authors. Using 
the WKB approach in potential theory, Bethe M derived, for proton-nucleus 
scattering, the expression 

$ = 21n(1.06/|ki|&e), (4) 

where |k x | is the cm. momentum, b is the range of strong-interaction forces 
defined by the size of the nucleus, and the cm. scattering angle. Similar 
results were derived within potential theory by several authors |J, ||. 
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(2) 



A relativistic derivation of the phase was attempted by Soloviev M , who 
obtained 

$ = 21n(2/0), (5) 

a result which differs considerably from the result of Bethe @] . Utilizing the 
same technique, Rix and Thaler [g] derived a result quantitevely close to that 
of Bethe [|. 

West and Yennie M , as already mentioned, obtained the phase of Coulomb- 
nuclear interference using Feynman diagrams. For a conventional parametriza- 
tion T N (s,t) ~ exp(— B\t\/2), the result of West and Yennie reads 



$>w-Y = T 



\n(B\t\/2)+ 1 + 0(B\t\) , (6) 



where 7 = 0.577... is Euler's constant. The upper (lower) sign corresponds 
to the scattering of pp (pp). 

Cahn || analized also the effect of the electromagnetic form factor and 
obtained a general expression for the phase. The results of |§ were in com- 
plete agreement with and especially the formula ([]) was derived from a 
rather different perspective. 

The main difference from the result M is a shift of the Coulomb ampli- 
tude due to the form factor's influence on the phase. If we introduce the 
electromagnetic form factor of the proton, the Born term of the Coulomb 
amplitude has the following form 

T C (s,t)=T 8 -^n\t\h (7) 

where the form factor may be chosen as 

f(\t\) = e ^ /A \ A 2 = 0.71GeV 2 . (8) 
In this case the Coulomb phase has the following form 

'B\t\ 



&Cahn — T 



In 



(^)+7 + ln(l + 



BA 2 , 

+ (4|t|/A 2 )ln(4|t|/A 2 ) + 2|t|/A 2 ]. (9) 

All these results were obtained under the assumption that \t\ — > 0. The 
derivation of the phase in a large domain of momentum transfer was at- 
tempted by Selyugin pj, and Kopeliovich and Tarasov |12|]. In the region 
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of interest for the present paper < \t\ < 0.01 GeV 2 , the latter results are 
similar to that of Cahn ||] and the main difference is in the dip region of the 
differential cross-sections. The phase obtained in [[| accurately takes into 
account the dipole form factor 



A 2 



LA 2 + \t\ 



A 2 = 0.71 GeV 2 



(10) 



and the complicated structure of the nucleon amplitude. 
The phase has the following form 



Selyugin nr 



In H + 2 7 
4 



57TS 



T N (s,t) 



bdbJ (bV^i)5 c (s } b) 



,2iS N (s,b) 



where 



0.111og(l + 400|t|) 



(11) 



takes into account influence of the form factor on the pure coulombic phase, 
and 



(12) 

The prescription of West and Yennie H] was successfully used by several 



6 c (s,b) = \nb + K (bA) + ^bkK^bk) 



authors [1C] and [] L| for describing differential cross sections in low t region. 

In what follows we will investigate four different cases of the Coulomb 
phase - the phase calculated with the nucleon amplitude of the model [|TJ 
(which does not acquire any additional parameter) with the prescription of 
West and Yennie 0, the phase calculated with prescription of Cahn ||, the 
prescription of Selyugin M, and the phase equal to zero. 



2 THE NUCLEAR AMPLITUDE 

We believe that any nuclear amplitude that is capable of a high accuracy 
description of the combined set of high energy pp and pp data (total and 
differential cross sections, p parameter etc.) over the entire \t\ spectrum, if 
properly combined with the correct Coulomb amplitude must account well 
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for the data in the interference region. That this is so, we will prove using the 
particular nuclear amplitude which has been derived in jlj to describe total 
and differential cross sections at high energies (y/s > lOGeV) in the range 
of momentum transfer 0.01 < \t\ < 14.5 GeV 2 using the eikonal approach 
(another one could have been the amplitude of fHI). We just write the 
nuclear amplitude of M 



T(s,b) 



3 2i5(s,&) 



1 



2i 



(13) 



where the eikonal has the following form 

5%{s, b) = 5+ (a, b) + 5+(s, b) + (J+ (s, b) T5®(s,b) + 5+(s, b) T 5~(s, b). (14) 

We refer the reader to the original literature for details; let us simply recall 
that here (s, b) are Pomeron contributions. The superscript '+' de- 
notes C even trajectories (the Pomeron trajectories have quantum numbers 
+ J ++ ), while ' — ' denotes C odd trajectories. Sq(s, b) is the Odderon con- 
tribution (i.e. the C odd partner of the Pomeron whose quantum numbers 
are _ J ); <5j and 8~(s,b) are the contributions of secondary Reggeons, (/ 
as representative of the C = +1 families and uo of the C = —1). 

In order to relate t- and 6-spaces one proceeds via Fourier-Bessel trans- 
forms 

poo 

db 2 j (bV^i)f(b) , 



fit) = 4ns 



(15) 



f(b) 



167TS 



dtJ (bV-t)f(t) 



Using this parametrisation we obtain the following expressions for the 
Coulomb phase 



£i*<(M = o) 



W-Y 



T- 



ln 



Mill* 

4 



) +7 



£i*M = 0) 



(16) 



and 



Cahn 



ln 



Pi(s)\t\ 
4 



+ 7 + In 1 + 



1(5 



E 4 *i(M=0) 

T (4|t|/A 2 )ln(4|t|/A 2 ) T 2|t|/A 2 



(17) 
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where p\ = 4c^(0) lns/so+rf, i = Pi, Pi, P3, O, f, uj. The upper (lower) signs 
correspond to pp (pp). 

Crossing symmetry is restored by replacing s — > se~ ln ^ 2 . We introduce 
the dimensionless variable 

s = -e-*5 , (18) 
so 

in terms of which we give each C+ and C— contribution with its appropriate 
signature factor in the form 



S + (s,b) 



5~(s,b) 
P 2 - 



so 4np 2 ' 

4a'(0)ln5 + r 2 , 

(C = +l); 

_ b ' 2 

__ i ? (0)-l^ 



(19) 



s Anp 2 
4a' (0)ln5 + r 2 , 
(C = -l). 



(20) 



3 RESULTS 

In [Q, the adjustable parameters have been fitted over a set of 982 pp and pp 
datafj of both forward observables (total cross-sections a tot , and p - ratios of 
real to imaginary part of the amplitude) in the range 8. < yfs < 1800. GeV 
and angular distributions (%) in the ranges 23. < ^fs < 1800. GeV, 0.01 < 
\t\ < 14. GeV 2 . A good x 2 /d.o.f. = 2.60 was obtained and the parameters 
are given in Table |l|. 



The data are available at 



REACTION DATA Database http://durpdg.dur.ac.uk/hepdata/reac.htmi 
CROSS SECTIONS PPDS database tyttp://wwwppdsAhep.su:8001/c5-5A.hlmi 



http:/ /pdg. lbl.gov/2000/ contents- plots. htm\ 
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Pomerorii 




f-Reggeon 




0.0578 ±0.0020 


A f 

J 


-0.31 (FIXED) 


QPi 


53.007 ±0.795 


c f 

J 


191.69 ±2.12 


aL 

Iri 

< TO) 


0.5596 ±0.0078 (GeV' 2 ) 
6.3096 ±0.2522 (GeV' 2 ) 


i 
r 2 f 


0.84 (GeV- 2 ) (FIXED) 
31.593 ± 1.099 (GeV~ 2 ) 




Pomerori2 




cu-Reggeon 




0.1669 ±0.0012 


A w 


-0.53 (FIXED) 




9.6762 ±0.1600 




-174.18 ± 2.72 


a™ 


0.2733 ±0.0056 (GeV- 2 ) 


a' 

Ul 


0.93 (GeV- 2 ) (FIXED) 
7.467 ± 1.083 (GeV' 2 ) 


^2 

/ TIE 

Jr2 


3.1097 ±0.1817 (GeV~ 2 ) 


r 2 , 




Pomeron 3 






A P3 


0.2032 ±0.0041 


so 


1.0 (GeV 2 ) (FIXED) 


Cp 3 


1.6654 ±0.0669 






«p 3 


0.0937 ±0.0029 (GeV- 2 ) 








2.4771 ± 0.0964 (GeV- 2 ) 








Odderon 






A 


0.19200 ±0.0025 






Co 


0.0166 ±0.0022 






«o 


0.048 ± 0.0027 (GeV- 2 ) 






r 2 


0.1398 ±0.0570 (GeV~ 2 ) 







Table 1: Parameters obtained in |l[]. 

We now consider the complete set of data including the Coulomb region 
which consists of 2158 dataQ. To these data, we apply the model including 
the Coulomb part with its phase and we simply plot the physical quantities 
without any additional fitting. 

The total cross sections and the ratios of real to imaginary parts of the 
forward amplitude are presented in Figs. [I] and |^. 

5 The data are available at 
REACTION DATA Database http://durpdg.dur.ac.uk/hepdata/reac.htm, 



CROSS SECTIONS PPDS database http://wwwppds.ihep. su:8001/c5-5A.htmi 



http:/ /pdg. lbl.gov/2000/ contents- plots. htmi 
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Figure 1: Total cross sections of 
pp scattering (hollow circles) and 
pp scattering (full circles) and cor- 
responding curves in the model of 
Ref.ffl. 




is (GeV) 

Figure 2: Ratios of the real to the 
imaginary part of the forward pp 
scattering amplitude (hollow circles) 
and pp scattering amplitude (full cir- 
cles) and corresponding curves in the 
model of Ref. 0. 



In order to compare different approaches to the Coulomb phase, we have 
calculated the \ 2 f° r the region of low \t\ : < \t\ < 0.01 (GeV 2 ) in four 
different cases: 

1. The Coulomb phase is equal to zero. 

2. The Coulomb phase is calculated with the prescription of West and 
Yennie ([161). 



3. The Coulomb phase is calculated with the prescription of Cahn (|17|). 

4. The Coulomb phase is calculated with the prescription of Selyugin ([TT|) 
The results may be found in Table |2|. 
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Coulomb phase 


Number of points 


X 2 per point 


$ = 


604 


3.49 


$W-Y (H§ 


604 


2.09 




604 


1.86 




604 


1.84 



Table 2: \ 2 P er point for the region of low \t\ : < |i| < 0.01 (GeV 2 ). 

As is seen from the Table |2|, the experimental data marginally "prefer" 
the Coulomb phase calculated with the prescription of Cahn |§ and Selyugin 
[§] over that of West and Yennie but taking the Coulomb phase equal to 
zero is excluded by the data and this is gratifying on physical grounds. The 
difference between the phases of [H] and |9j is negligible in the small \t\ region 
of the interest (as may be seen in Fig. [jD, though the exact result obtained 
in |J over the whole kinematical region of \t\ shows that the difference is 
drastic in the region of the diffractive dip (as seen in Fig. ||). 

Apart from this, the Coulomb phase calculated in the Selyugin approach 
exhibits a non trivial behaviour of the real part (which has a zero in the point 
of the diffractive dip), and of the imaginary part (which has a minimun in 
the same point) Fig. f|. 

The diffrence between the approaches of Cahn || and of West and Yennie 
[0] is more pronounced even in the region of low t as may be seen in Fig. [3[ 

It is of some interest to perform a more detailed comparison of the \ 2 
derived for the region of low momemtun transfer. This can be found in 
Table | for $ W -y, ®Cahn, and $ S eiyugin- 

With all due caution, a general pattern emerges; some data points have 
anomalously large \ 2 values and these are the same in all cases. 

The pp differential cross section at y/s = 546 GeV, for instance, sticks 
out conspicuously. The explanation is the normalisation of the cross section. 
The sysytematical error allows to change the normalisation of the data in 
the ±10% corridor. If we could multiply our predictions by a factor 1.06, we 
would have a good description of this set of data, as is seen in Fig. || and |[ 
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52.8 GeV 



52.8 GeV 





-t (GeV ) 



Figure 3: The Coulomb phase calculated 
in the framework of Selyugin |9j (the solid 
line) and of Cahn (the dotted line) and 
of West and Yennie Q (the dashed line) in 
the region of small t. ^/s = 52.8 GeV 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 

n . . -t(GeV 2 ) 

r lgure 4: The Coulomb phase calculated 
in the framework of Selyugin j^] , real part 
(the solid line) and imaginary part (the 
dashed line) and in that of Cahn Q| (the 
dotted line) in the region of high t. 
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r lgure 5: Differential cross section of pp 
scattering and curves corresponding to its 
description in the model. The dashed line 
represents our predictions. The solid curve 
is our predictions normalized by a factor 
1.06. 
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Figure 6: The same as in Fig. || for low 
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Data 




< \t\ < 0.01 (GeV 2 ) 
#points, ntot 


XVntot 

&W-Y 


xVntot 

<&Cahn 


X^/ntot 

Selyugin 


1 

2 
3 
4 


pp 

" total 

PP 

total 
p PP 

p p P 




33 
68 
11 
48 


0.7167 
0.3617 
0.6086 
0.4326 


0.7167 
0.3617 
0.6086 
0.4326 


0.7167 
0.3617 
0.6086 
0.4326 


5 
6 
7 
8 
9 

10 


daPP 
dt 

daPP 
dt 

do-PP 
dt 

do-PP 
dt 

daPP 
dt 

do-PP 
dt 


24.3 

30.4 

52.6 

62.3 

541.0 

546.0 


26 
21 
12 
3 

37 
15 


1.4965 
1.4702 
1.9847 

I. 0263 
2.0811 

II. 4692 


2.4879 
1.4130 

I. 4676 
0.8135 
2.3993 

II. 0364 


2.8105 
1.4027 
1.3440 
0.7635 
2.5246 
10.9008 


1 1 

± L 

12 
13 
14 
15 
16 
17 
18 
19 
20 
21 


do-PP 
dt 

doPP 

toPP 

toPP 
dt 

doPP 
dt 

doPP 
dt 

doPP 
dt 

doPP 
dt 

doPP 
dt 

doPP 
dt 

doPP 
dt 


lu.U 

12.3 
19.4 
22.2 
23.9 
24.3 
27.4 
30.7 
44.7 
52.8 
62.3 


oo 

37 
45 
45 
94 
25 
40 
8 

25 
21 
4 


1.5731 
1.9674 
1.4666 
2.3910 
1.1418 
2.2246 
0.3174 
5.0147 
7.1074 
8.5177 


9 D9Q7 
1.4169 
1.7840 
1.3374 
2.3699 
1.2638 
2.0014 
1.2071 
2.9964 
4.2590 
6.7054 


1 qri n 

-L . yuiu 

1.3931 
1.7509 
1.3192 
2.3087 
1.3587 
1.9607 
1.5482 
2.7679 
3.9018 
6.2988 



Table 3: y 2 per point for the region of low \t\ : < \t\ < 0.01 (GeV 2 ) for 

&W-Y, &Cahn, and $Selyugin- 



To appreciate the role of the Coulomb phase, we plot the angular distri- 
butions with the appropriate Coulomb phase and with the phase taken equal 
to zero in the following figures [7], |^. Even though the difference is minimal, 
the numerical conclusion is that the data, quite unambiguously, prefer the 
appropriate nonzero Coulomb phase. 

The comparison with the West-Yennie and Cahn phases is not shown 
because the plot would be undistinguishable from the Selyugin one. 

Let us remind that no additional fitting was done. 

We do not report here any of the pp and pp angular distributions which 
one obtains over the full range of \t\ values because the original reproduction 
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of these quantities is left basically unchanged by the Coulomb amplitude and 
the interested reader is referred to the original paper [|l|. We simply report 
here, for completeness, the predictions at RHIC and LHC energies both in 
the interference region and over the entire \t\ range. Predictions of the model 
and comparison with the nuclear amplitude for RHIC and LHC are shown 
in figures f 



10, 11, 12 



11)01) 

'too 
soo 

700 
, 600 
500 



• 
























I 












I 






1 






\ 


22.2 GeV 




\ 

\ 





































B 700 

iS 600 



0.015 

-t (GeV 2 ) 





















































30.4 GeY 




X 

K 

























0.015 

-t (GeV 2 ) 



Figure 7: Differential cross section 
of pp scattering and curves corre- 
sponding to its description in the 
model. The solid line corresponds to 
the Coulomb phase calculated with 
prescription of Selyugin flTID and the 
dashed line corresponds to the zero 
phase. 



Figure 8: Differential cross section 
of pp scattering and curves corre- 
sponding to its description in the 
model. The solid line corresponds to 
the Coulomb phase calculated with 
prescription of Selyugin (|TT|) and the 
dashed line corresponds to the zero 
phase. 
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2 4 6 8 10 12 14 



-t (GeV 2 ) 

Figure 9: Prediction of the model for the pp scattering differential cross 
section at RHIC (y/s = 100, 500 GeV) and LHC (s/s = 14 TeV). A 10" 4 
factor between each successive set of data is omitted. The prediction for 
RHIC at the energy of y/s = 500 GeV is multiplied by a factor of 10~ 4 and 
the prediction for LHC at the energy of y/s = 14 TeV is multiplied by a 
factor of 10" 8 . 
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Figure 10: Prediction of the model for the pp scattering differential cross sec- 
tion at RHIC (s/s = 100 GeV). The solid line corresponds to full amplitude 
including the Coulombic one and the dashed line corresponds to the nuclear 
amplitude M. 
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Figure 11: Prediction of the model for the pp scattering differential cross sec- 
tion at RHIC (y/s = 500 GeV). The solid line corresponds to full amplitude 
including the Coulombic one and the dashed line corresponds to the nuclear 
amplitude M. 
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Figure 12: Prediction of the model for the pp scattering differential cross 
section at LHC (s/s = 14 TeV). The solid line corresponds to full amplitude 
including Coulombic one and the dashed line corresponds to the nuclear 
amplitude M. 



To show the influence of the interference between nucleon and Coulomb 
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amplitudes we show the ratio 



R 



dt dt 



d,cr r 



dt 



da_ 

dt 



(21) 



in figures and [L4| and we conlude that the interference term becomes 

negligible at RHIC and LHC energies at t 0.01 GeV 2 for RHIC and 

t 0.005 GeV 2 for LHC. 
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Figure 13: The interference term 
between hadron and Coulomb am- 
plitudes for at the RHIC energy. 
The dashed line corresponds to the 
Coulomb phase calculated with the 
prescription of Selyugin and the 
dotted line to the zero Coulomb 
phase. 
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Figure 14: The interference term 
between hadron and Coulomb am- 
plitudes at the LHC energy. The 
dashed line corresponds to the 
Coulomb phase calculated with the 
prescription of Selyugin |J and the 
dotted line to the zero Coulomb 
phase. 



Now let us return to the problem of extracting the ratio p(s, t = 0) of the 
real to the imaginary parts of the scattering amplitude. A general theorem 
[ |i"4jl shows that this ratio as a function of the momentum transfer could not 
have a constant sign in a strip sm < s < oo, — T < t < 0, for any s M and 
T > 0. The method of extracting the ratio includes a simplification, the 
scattering amplitude is usually presented in the following form: 
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T(s, t) = T ^f(\t\) + (* + p(s, 0))sa tot e Bt (l - za$) (22) 

where the ratio p{s,t) is approximated by its value at t — 0. Should this 
simplification affect the measurement? The answer would be "yes" should 
the change of sign of p occurr close to the region of the maximum of the 
Coulomb nuclear interference. 

In our particular model, the zero of the real part of the scattering am- 
plitude goes to zero with the increase of energy, but remains far from the 
region of the maximum of the interference. For example, at LHC energies 
^fs = 14 TeV, the maximum occurs at \t max \ ~ 10~ 3 GeV 2 (fig.|i~4"D while the 
zero is at \t zero \ ~ 0.16 GeV 2 ( fig. ITS], 116]), so the situation is not too bad. 



1U 1U 1U 1U 

Vs (GeV) 

Figure 15: The position of the 
zero of the real part of the scatter- 
ing amplitude in a three-component 
Pomeron model |T[| (solid line) and 
in a "geometrical scaling" model 
tzero = -l/(Alns), A = 0.25 GeV~ 2 
[JLJ (dashed line). 
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Figure 16: The real part of the scat- 
tering amplitude at the LHC energy 
in the model (Jj. 
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CONCLUSION 



All the three choices of the Coulomb phase give good description of the 
existing data, though in terms of \ 2 P er point the data "prefer" the phase 
calculated with the prescriptions of Cahn || (|I7D and Selyugin || (fLip . We 



conclude that such prescriptions are a good basis for the evaluation of the 
Coulombic contribution in the full scattering amplitude. We expect that 



similar conclusion would apply to an analysis repeated along the lines of |T2 
but we have not performed this analysis. 

As we have seen, the addition of the nuclear amplitude (with parameters 
fitted from total and differential cross sections) and of the Coulomb one 
(with its proper phase) gives a total amplitude which reproduces quite well 
the data in the interference region without any additional parameters and 
with no need to refit existing ones. 

This allows us to predict the RHIC Coulomb interference which requires 
the measurements to start from \t\ < 0.005 GeV 2 at the energy of ^/s = 
100 GeV and from \t\ < 0.004 GeV 2 at the energy of y/s = 500 GeV. 
Likewise, LHC will be able to cover the Coulomb region if the measurement 
starts from \t\ < 0.001 GeV 2 . 
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